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Draw A Sphere (input Radius, CenterX, CenterY, VectLight, VectViewer, 

VectorPole, VectorEquator) 
// the image inputs include the size of the sphere, where it is to be drawn, 
// where a lighting source is positioned and where a viewer is positioned 

{ 

// set up initial vectors 

vectSpecularHighlight = Normalize(vectViewer + vectLight); 
vectPoleCrossEquator = VectorPole cross VectorEquator; 

// prepare lookup tables, can be computed before rendering 

// portions of later calculations pre-calculated here b/c x & y invariant to 

// other parameters 

for ( i = -rad; i <= rad; i++) 

{ j = i * 1 / rad; 

xMultiplyDiffuse[i] = j * vectLight.x;; // setup diffuse component 

yMultiplyDiffuse[i] = j * vectLight.y; 

xMultiplySpecular[i] = j * vectSpecularHighlight.x; // setup specular 
yMultiplySpecularji] = j * vectSpecularHighlight.y; 
xMultiplyPole__LUT[i] = j * vectorPole.x; // used for texture 
yMultiplyPole_LUT[i] = j * vectorPole.y; 
xMultiplyEquator_LUT[i] = j * vectorEquator.x; // setup equator 
yMultiplyEquator_LUT[i] = j* vectorequator.y; 
xMultiplyPE_LUT[i] = j * vectPoleCrossEquator.x; // where pole & 
xMultiplyPE_LUT[i] = j * vectPoleCrossEquator.y; //equator cross 

} 

for ( x = 0; x < rad; x++ ) // finite set of discriminants 
{ disc = r A 2 - x A 2; 

for ( y = 0; y < x; y++ ) // thus finite set of z values 

{ disc = disc - y A 2; 
if ( disc > 0 ) 

{ zlnvRad = 1 / (squareroot(disc) * radius; 
zMultiplyDiffuse_LUT[disc] = zlnvRad * vectLight.z; 
zMultiplySpecular_LUT[disc] = zlnvRad * 

vectSpecularHighlight.z; 
zMultiplyPole_LUT[disc] = zlnvRad * vectorPole.z; 
zMultiplyEquator_LUT[disc] = zlnvRad * vectorEquator.z; 
zMultiplyPE_LUT[disc] = zlnvRad * 
vectPoleCrossEquator.z; 
} // end if 
} // end for y 
} // end for x 
// proc confd on Fig. 7b 



FIG 7A 



// Iterate over the scanlines in the sphere 
// combining the precomputed lookup elements as you go 
// for each scan line 
for ( y = -rad; y <= rad; y++ ) 
{ RadiusSubYSquare = r A 2 - y A 2; 

Bound = edgeBuffer[abs(y)]; // bound is the horizontal displacement from 

// y axis 

for ( x = (-bound + 1); x <= bound; x++ ) 
{ 

152 // iterate over every pixel in the scanline y 

disc = RadiusSubYSquare - x A 2; // compute disc for look up table 

// index 

diffuse = yMultiplyDiffuse[y] + xMultiplyDiffuse[x] + 

zMultiplyDiffuse_LUT[disc]; 
specular = yMultiplySpecular[y] + xMultiplySpecular[x] + 

zMultiplyDiffuse_LUT[disc]; 
specular = SpecularRemapLUT[specular]; // remap to range 0 -1 .0 

150 . 

// compute the final intensity for a pixel 
intensity = diffuse * diffuseFactor + specular * specularFactor + 

ambient * ambientFactor; 
// compute the u & v texture components for a pixel 
NormalDotPole = xMultiplyPole_LUT[x] + yMultiplyPole_LUT[y] + 

zMultiplyPole_LUT[z]; 
NormalDotEquator = xMultiplyEquator_LUT[x] + 

yMultiplyEquator_LUT[y] + zMultiplyEquator_LUT[z]; 
latitude = arccos(NormalDotPole); 
vTexture = latitutde/PI; 

longitude' = NormalDotEquator / sine(latitutde); 
clamp longitude 1 to range -1.0 to 1.0 
longitude = arccos(longitude'); 

// determine how longitude wraps around sphere 
if (xMultiplyPE_LUT[x] + yMultiplyPE_LUT[y] + 

zMultiplyPE__LUT[disc] < 0 ) 
{ uTexture = longitude; } 
else 

{ uTexture = 1 - longitude; } 

//fetch a textured pixel from coordiante uTexture, vTexture 
// scale intensity of textured pixel by Intensity 
// draw the lit, texture pixel at location ( x + centerX, y + centerY) 
} // end for x 
} // end for y 

}//endproc FIG 7B 
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